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ABSTRACT 

The optical spectra of objects classified as QSOs in the SDSS DR6 are an- 
alyzed with the aim of determining the value of the fine structure constant in 
the past and then check for possible changes in the constant over cosmological 
timescales. The analysis is done by measuring the position of the fine structure 
lines of the [OIII] doublet (AA4959 and AA5008) in QSO nebular emission. From 
the sample of QSOs at redshifts z < 0.8 a sub sample was selected on the ba- 
sis of the amplitude and width of the [OIII] lines. Two different method were 
used to determine the position of the lines of the [OIII] doublet, both giving 
similar results. Using a clean sample containing 1568 of such spectra, a value of 
Aa/a = (+2.4 ± 2.5) x 10~^ (in the range of redshifts z ~ — 0.8) was deter- 
mined. The use of a larger number of spectra allows a factor ~ 5 improvement 
on previous constraints based on the same method. On the whole, we find no 
evidence of changes in a on such cosmological timescales. The mean variation 
compatible with our results is 1/ < t > Aa/a = (+0.7 ±0.7) x 10"^'' yr"^- The 
analysis was extended to the [Nelll] and [SII] doublets, although their usefulness 
is limited due to the fact that all these doublets in QSOs tend to be fainter than 
[OIII], and that some of them are affected by systematics. 

Subject headings: atomic data-line identification-line profiles-quasars: emission 
lines 



1. Introduction 



In Kaluza-Klein and super string theories that try to unify all the fundamental inter- 
actions, the constants of nature are functions of a low mass dynamical scalar field, which 
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slowly changes over cosmological timescales ( lUzan 1 12003| ). Astrophysics offers a possible 
test to constrain the parameters of such theories by directly measuring the values of such 
constants through the comparison of properties of objects at different evolutionary epochs of 



the U ni verse. Reviews of the techniques used have been presented by iLandau and Simeone 



( l2008l ). iKanekar et al. I (120081 ) and others. The biggest observational effort to constrain the 
values of physical constants in the past has been directed at the fine structure constant 
(a = e'^/hc), which plays a fundamental role in the characterization of electromagnetic in- 
teraction. Several methods have been dev eloped based on the analys is of cosmic microwav e 
background data (INakashima et al. II2008I ). Big Bang nucleosynthesis ( llchikawa et al. II2002I ). 
and the fine structure splitting of several atomic lines in QSO sp ectra (see below). A review 
on the techniques and observational constraints can be found in iGarcia-Berro et al.l (120071 ). 



Here, we mention only those methods based on fine structure splitting. The splitting 
ratio (A2 — Ai)/(A2 + Ai) at t wo different epochs gives the relative difference in a between 
these two epochs. It is shown (lUzan II2003I ) that 



Aa 



[(A2-Ai)/(A2 + Ai)], 
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a ■ 2 l[(A2 - Ai)/(A2 + Ai)]o 
where A2 and Ai are the wavelengths of the pairs of the doublet, and the subscripts z and 
refer to the values at redshift z and local respectively. 

The method has been employed to measure the r elative separation of absorption lines 
in the spectra of QSOs, inc luding the alkali do ublet ( iBahcall et al. 1119671 ). and the many 



multiplet method (MML), (IWebb et al. Ill999[). The best constraints obtained using the 



alkah doublet method are those by lChand et al. I (120051 ) [Aa/a = (0.15 ± 0.44) x 10 ^] over 



the range 1.59 < z < 2.32]. The MML simultaneously analyzes many doublets of many 
atomic species, and due to the large number of lines used, higher precision than with the 
alkali doublet method. 

T he MML is the only method that has res ulted in claims for the detection of variation 
of a JWebb et al. I Il999l : iMurphv et a"l~l I2OO3I ) Aa/a = (-0.5 74 ± 0.102) x 10"^ in the 



range 0. 2 < z < 3.7, a . 



see also (IMurphv et al. 



thoug h these results are controversial ( IChand et al. 1 120061 ) , (but 



20081) on the opposite view), and are at odds with similar studies 



dChand et aLll2004i ISrianad et a"l~l 120041 : Ihevshakov et al. II2OO7I ) 



The MML method suffers from several uncertainties a nd might be affected by possible 

( 120041 ) among others. Another group 



systematics, as was poi nted out by IBahcall et al. 



( iLevshakov et al. Il2005[ ) has developed a slight modification of the MML method in which 



only one atomic ion (Fe II) is used, avoiding many of the assumptions and uncertainties 
inherent to the MML method. Their latest analysis (ILevshakov et al. 1120071 : iMolaro et al. 
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20081) determined Aa/a = (-0.12 ±1.79) x lO"*' at z = 1.15 and Aa/a = (5.4 ±2.5) x 10"^ 
at z = 1.84 although the authors cautiously stated that there could be some unnoticed 
systematic effects that might challenge that interpretation. 

In this paper, we use the [OIII] neb ular emission lines in QSOs t o constrain past varia- 
tio ns in a. The method w as p roposed bvlBahcall and Salpeter I (119651 ) , and was later applied 
by iBahcall et al. I (120041 ) and iGrupe et al. I ( 120051 ) . Because the analysis is based on a pair 
of lines only, the constraints on Aa/a are not as strong as those obtained with the other 
methods mentioned above, but has the advantage that it is more transparent and less subject 
to systematics. The method is based on the same element and level of ionization, and both 
lines originate in the same upper energy level, so the analysis is quite independent of the 
physical conditions of the gas where the [OIII] lines originat e. It therefore r e prese nts a good 



alternative for determining the value of a on a firm basis. iBahcall et al. I (120041 ) analyzed 



the QSOs of the Early Data Release of SDSS and built a clean sample of 42 QSOs in the 
range 0.16 < z < 0.80, from which they derived /S.a/a = (±0.7 ± 1.4) x 10^'^. Here, we use 
the QSOs in the latest release of the SDSS to improve such constraints significantly. 



2. Sample selection and methodology 

The SDSS-DR6 Catalog Archive S'erwe?]^ contains 77082 objects classified as QSOs. By 
SQL queries we downloaded the extracted Id spectra of all of them. These spectra have 
a resolution ~ 2000 and the range covered is ^3400-9200 A, which includes the [OIII] 
doublet! up to redshift ^ 0.85; this restriction in redshift automatically limits the sample 
to 28,860 spectra. The SDSS pipeline also provides continuum-subtracted spectra, which 
in most cases is quite acceptable; however, we found some cases (usually related with the 
presence of especially strong and wide Balmer lines) that the estimation of the continuum 
was unsatisfactory. To proceed in a consistent way we decided to recalculate the continuum 
by fitting a cubic spline to all the 1-d spectra, masking those regions that were centered 
around the nominal rest frame wavelengths of strong emission lines. We inspected several 
hundred spectra chosen randomly, and checked that this new estimation of the continuum was 
reasonably good. Anyway, we have checked that the restrictions found for Aa/a (see next 
sections) do not change significantly if we use, instead of our estimation of the continuum, the 
one provided by SDSS. The wavelength calibration was checked by measuring the position of 
the atmospheric 01 line. From 1656 sky spectra we determined a mean position of 5578.885 



^http: / /www. sdss.org/ dr6 / access /index. html#C AS 
^ Other species are considered in Section 3. 
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A and a standard deviation of 0.287 A which agrees quite well with the theoretical value 
(5578.887 A) and with expectations from the spectral resolution of the spectra. The selection 
of the final sample was done on the basis of the strength, shape and width of the [OIII] 
lines, reliability of the continuum estimation, and possible contribution due to the relative 
proximity of the H/3 line. The estimation of the spectral position of the [OIII] lines was done 
following two methods which are described below. 



2.1. Method 1 



To estimate the centroids of the lines, we first determined the FWHM of each line 
and averaged the flux of the pixels within the spectral range covered by that FWHM. We 
estimate the relative strength of the [OIII] lines with respect to the noise in the adjacent 
continuum. Such noise was estimated by measuring the continuum rms of each spectrum in 
the wavelength range 5040-5100 A (rest frame). We select those objects that have strong 
{SNR > 20 in the peak) [OIII] lines. These restrictions provide us with a sub sample (which 
we refer to hereafter as the "raw sample") of 3739 objects. The contribution of the H/3 
line could in principle produce a blueshift in the estimation of the centroids of the lines, 
particularly affecting the 4959 A line. We checked that this contribution could be properly 
quantified by computing the level of the flux (after subtraction of the continuum) in a spectral 
region slightly bluer than the position of the [OIII] (AA4959) A line. After many trials we 
chose the 4925-4935 A (rest frame) wavelength range to estimate possible residuals of H/3 
and remove from the sample those objects with fluxes in that interval above 0.05 the peak 
value of the [OIII] (AA4959) A line. Fig{T] shows an example of a spectrum that has been 
rejected on the basis of that criteria. We also eliminated spectra with very wide or double 
peaked [OIII] lines. Figure [2] shows the different spectral regions used to calculate fluxes 
and centroids, and to assess potential HP residual contamination and the continuum level. 
We do not impose any further restrictions based on the shape (the presence of asymmetries) 
of the [OIII] lines. After all these cuts, 1978 spectra remained to constitute our "clean 1" 
sample. Using this sample we obtain mean Aa/a = (—0.9 ± 2.6) x 10~^. The results are 
very robust against the precise constraints; for instance selecting only those spectra with 
SNR > 30 in the peak of the [OIII] lines, we obtain mean Aa/a = (1.9 ± 2.6) x 10"^ 
Figure [3] shows examples of the region centered around the [OIII] doublet of some spectra 
randomly chosen from this sample. The figure shows that for a given spectrum both lines 
are well above the noise level, and have similar shapes. A few of them show the presence 
of asymmetries probably related to the kinematics of the cloud or the presence of multiple 
clouds. The distributions of centroids of each member of the doublet show similar rms (0.78 
A); this indicates the absence of significant distortions in the estimation of the position of the 
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4959 A line due to the relative proximity of the H/3 line. The distribution of the wavelength 
separation between the centroids of both lines has a rms of 0.11 A; this indicates that the 
main factor broadening the distribution of the centroids of the [OIII] lines are absolute 
spectral shifts which mostly reflect the kinematics (< 150 km s~^) of the narrow line clouds 
where the [OIII] originates. 

2.2. Method 2 

The emission lines of OIII (4959 and 5008), and H/3 (narrow and broad components) 
were simultaneously modeled by single Gaussian each. The wavelength position of each line 
of the [OIII] doublet was directly estimated as the central position of the corresponding 
Gaussian. In principle this seems a very simple description of the line profiles, and in many 
cases produces a poor fit to the data. However, in practice the method takes advantage of the 
expected similar shape for both lines of the doublet, and removes most of the contribution 
of H/3 in the spectral region of the [OIII] doublet. After removing the most extreme cases 
of poor fits, the results are quite consistent and robust. Selecting those spectra in which 
each of the [OIII] lines are described by Gaussian with peak amplitudes > 20 the level of the 
continuum noise, and a within the range 1.4-3.0 A (rest frame of the [OIII] 5008 line) there 
were 1568 spectra ('clean sample 2') from which we obtained Aa/a = (+2.4 ± 2.5) x 10~^. 
Figure H] shows the wavelengths determined for the spectral position of both lines of [OIII] 
(to compute the rest frame wavelengths we have used the estimates of redshifts provided by 
the SDSS database). The wavelengths of both lines lie along a line crossing from the bottom 
left to the top right and don't show any systematic effect. 

3. Results and discussion 

The results obtained by both methods are compatible and quite similar. Both estima- 
tions are also compatible with the local value. Most of the spectra included in the clean 
sample 2 (1223 out of 1568) are also included in the clean sample 1. 

Fig O shows a comparison between the results of both methods. There is a general good 
agreement between them, although the method of fitting Gaussian tends to be a bit more 
accurate (the rms of the distributions of Aa/a are 9.8 x 10~^ and 8.0 x 10~^ for the clean 
samples 1 and 2 respectively). Then, here after along the paper, we will use the results 
obtained by the second method. 

To analyze the value of Aa/a as a function of redshift (or look-back time) we compute 
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the mean values in redshift interval. This is shown in Fig. |6] and Table 1, which are the main 
results of this paper. The bins in redshift have been built in order to include approximately 
the same number (175) of spectra each. Although none of the bins shows a statistical 
significant departure from zero, the largest deviation is at the 2.1a level in the bin which 
corresponds to the range in redshift 0.110-0.152. That bin includes those cases in which one 
of the [OIII] lines lies near the spectral position of the atmospheric 01 line. To check for the 
possible influence of some residual of the 01 line in the estimation of the centroids of [OIII] 
in that range of redshift, we built a new sample excluding those spectra having redshifts in 
which any of the [OIII] lines is closer than 15 A to the spectral position of the atmospheric 
01 line. This new constraint removes 59 spectra. The overall results on Aa/a do not change 
much (Aa/a = (+0.6 ± 2.4) x 10~^) using this new clean sample, but the significance of 
the departure from zero in the corresponding bin in redshift (0.107-0.161) is largely reduced 
(Aa/a = (+7.0 ± 5.8) x 10"^). 

The range in redshift spanned by our sample corresponds to a maximum look-back time 
of 6.6 Gyr, and a mean of 3.0 Gyr. Following Bahcall et al., the mean rate of possible 
changes is 1/ < t > Aa/a = (+0.7 ± 0.7) x 10~^^ yr"^- A linear fit of Aa/a with respect 
to look-back time gives a slope of (0.2 ± 1.6) x 10^^^ yr~^- All these numbers as a whole do 
not show any significant evidence of change in Aa/a with redshift. 



3.1. Line ratios 

As a by-product of the analysis we have estimated the value of the line ratio of [OIII]: 

A(5008)AE(5008) 



A cx 



A(4959)AE(4959) 



where A and AE are the transition rates and the energy differences, respectively. The 
determination of this ratio is interesting as an observational test of calculation s from atomic 



theor y where a ~ 3 % level of discrepancy exists among the published results ( iGrupe et al 



20051 ). We performed numerous tests checking the dependence of A on possible residual 
contamination from HP (see previous section), SNR and [OIII] width. Our best estimation 
is obtained with the same constraints as those used for the estimation of Aa / a using method 
2. Fig. [7| shows the histogram of A values for the 1568 spectra included in 'clean sample 
2'; the resulting distribution is symmetric, centered at ~ 3 and has a standard deviation of 
0.10. From that distribution was obtained our best estimation of A = 3.012 ± 0.003 ± 0.010, 
where the error bars correspond to the la statistical and systematic errors respectively. The 
systematic errors correspond to uncertainties due to the criteria for selecting the samples. 
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The value obtained agrees w ith the estimate by Bahc all et al. (2.99 ± 0.02), with the 
measurements 2.993±0.014 bv lDimitriievic et aL I (120071) using 6 2 AGN, and it is compatible 
with theoretical estimates (2.98) by lstorey and Zeippen (2000). 



3.2. Other doublets 

As pointed out by Bahcall et al. and developed and used by Grupe et al., it is possible 
to constrain the value of a at cosmological look-back times using doublets of other emission 
lines that are common in the spectra of QSOs. The most promising are [NeV] (AA3346, 3426), 
[Ne III] (AA3869,3968), [01] (AA6300, 6365), and [SII] (AA6717, 6731). The main advantages 
of the use of these lines is the possibility of extending the study to higher redshifts and of 
checking the consistency of the results using the [OIII] lines. All of them lie in the optical 
range and, in the wavelength range of SDSS, can be useful for exploring the ranges in redshifts 
up to 1.69, 1.31, 0.44 and 0.37 respectively. However, these lines are usually weaker than 
[OIII] in QSO spectra and some of them also have drawbacks that limit their suitability for 
these studies. For instance, the lines of the [SII] doublet are relatively close so that further 
constraints on the width of such lines need to be applied. This pair is also relatively close 
to the usually very strong and wide Ha line in QSOs. On the other hand, the pair of lines 
for [NeV] and [Nelll] are separated by ~100 A and then spectral calibration uncertainties 
might be a more serious concern. 

We follow a similar procedure to the one used for [OIII], adapting the constraints on 
the SNR, possible contributions from neighboring lines, etc., to each specific doublet. As 
expected, because of the relatively low amplitude of such lines with respect to [OIII], the 
clean samples of each of these pairs are much smaller, and the sensitivity achieved is poorer 
than with [OIII]. We decided not to use the [NeV] and [01] doublets for this analysis because 
the number of spectra with lines strong enough for a meaningful analysis is very small. In 
order to have a representative sample of [Nelll] and [SII] we relax our selection criteria to 
include in the samples all the spectra having SNR> 10 in the peak of each line of the doublet. 
Table 2 presents the results for [Nelll] and [SII]. The values found from [SII] are compatible 
with those from [OIII]. Although the clean sample used in the case of the [SII] contains a 
relatively high number of spectra (466), the sensitivity achieved is comparatively much lower 
than the one obtained from [OIII] and [Nelll]. After doing exhaustive tests we found that 
the relative separation with respect to the width of the lines is the most important factor 
which limits the use of this pair. Figure |8] {left) shows a comparison between the Aa/a 
values obtained in common spectra between [SII] and [OIII] samples. The figure shows the 
compatibility of both estimations and the comparatively much large dispersion obtained 
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from the estimations based on [SII]. 

The values of iS^aja obtained from [Nelll] are significantly different from zero. We 
have checked that this departure from zero keeps along all the range in redshift analyzed 
using that doublet. The wide range in redshift covered by this sample and the absence of 
the effect for other lines (for instance, the lines of the [OIII] doublet or the atmospheric 01 
line) precludes possible systematic errors in wavelength calibration. Figure[8] (ri^f/it) shows a 
comparison between the estimations of Act/a from the [OIII] and [Nelll] respectively using 
only common spectra to both samples. Clearly, the results from [Nelll] are not compatible 
with those from [OIII]. For the distribution of the centroids of each member of the [Nelll] 
doublet, we have estimated mean values of 3869.61 and 3969.05 A and rms of 0.51 and 
0.58 A respectively. The comparison with the theoretical local values (3869.86 and 3968.59 
A respectively) shows discrepancies of 0.25 and 0.46 A to the blue and to the red respectively. 
If we restrict the sample according to stronger constraints on the SNR of both lines, these 
shifts tend to decrease only for the line at 3870 A. Although this could indicate inaccuracy in 
the tabulated wavelengths of both transitions, we think the main reason is that the estimated 
wavelengths for the line [Nelll] (3970 A) are affected by the proximity of several lines like He 
and the absorption line of Call-H. Nevertheless, if this is the case it would remain unclear 
why the rms of the the distributions of both members of the doublet are so similar, and why 
the shift in the 3970 A line does not decrease when only spectra with relatively high SNR 
are considered. Curiously enough, the high value found for Aa/a using [Nelll] is similar to 
the mean value found by Grupe et al. using this pair. 



3.3. Comparison with previous studies 

In Section 1 we have mentioned the different astronomical methods used to constrain 
the variation in a with particular emphasis on those based on fine structure splitting. Here, 
we make comparisons only with previous observations based on the relative splitting of pairs 
of emission lines in QSOs or AGN in general. The analysis by Bahcall et al. is based on 
the Early Data Release of SDSS, which comprises a sub sample of the sample analyzed in 
this paper, so it is worth making a comparison of both studies. There are 38(10) shared 
objects between the sample of these authors and our raw(clean) sample. Figure [10] shows 
a comparison between our estimates of Aa/a and those by Bahcall et al. The comparison 
shows a very good agreement between both studies. 

The distribution of Aa/a obtained from objects in common between our raw sample 
and the Bahcall et al. sample is similar, while the rms of Aa/a from our clean sample is 
0.5 X 10~^, slightly better than the sample of Bahcall et al. (0.7 x 10~^). So, although 
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the method and criteria for selecting the samples are quite different from each other, both 
studies agree quite well. The much higher number of spectra used in this study (1568 in our 
clean sample) with respect to the 42 objects in the Bahcall et al. sample would allow us in 
principle to reduce the uncertainty by a factor a/1568/42 = 6.1. The result presented by 
Bahcall et alE is (0.7 ± 1.4) x 10"^ , and therefore our result (Aa/a = (+2.4 ±2.5) x 10"^) 
improves the sensitivity by a factor ~5, which agrees with statistical expectations. 



Grupe et al. I (120051 ) applied the same method to a sample of 14 Sey 1.5 galaxies up to 
redshift 0.281. These authors include in their analysis the [Nelll], [NeV], [OIII], [01], and 
[SII] doublets. They found Aa/a = (1.5 ± 0.7) x 10~^, so the result presented in this study 
improves that sensitivity by a factor ~ 30, and the range in redshift covered by a factor ~ 3. 
The method used in this paper to constrain a in the past is about one order of magnitude 
less sensitive than methods based on the analysis of absorption systems along the line of 
sight of QSOs. However, as has been discussed above, the method is more transparent, does 
not require assumptions on the spatial distribution and kinematics of the gas involved, and 
in general is less subject to systematic errors. 



4. Conclusions 



The main findings of this work are: 



• We have analyzed the sample of 77092 objects classified as QSOs in the SDSS DR6, 
and from these we have selected a sample of 1568 objects with strong [OIII] emission 
lines up to redshift 0.80. 

• From the relative spectral position in such spectra of the nebular emission of [OIII] 
AA5007 and 4959 we have estimated a mean value for the fine structure constant of 
Aa/a = (+2.4 ± 2.5) x 10"^ through the last 6.6 Gyr. 

• Our analysis improves by a factor ~ 5 — 6 the existing uncertainty in the estimates in 
the possible change of a using the same method. 

• The maximum rate of change in Aa/a allowed by our analysis is (+0.7 ± 0.7) x 10~^^ 
yr~^ in the last 6.6 Gyr. 



•^In the appendix of the BahcaU et al. paper the SDSS Data Release One was analyzed and a value 
Aa/a (1.2 ± 0.7) x IQ-"* found. 
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Fig. 1. — Example of one spectrum that has not been included in the 'clean sample' on 
account of its broad H[5 line. 
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Fig. 2. — Diagram showing the regions used to compute the centroids and relative fluxes, and 
the regions used to assess the if/3 residual contamination and the noise in the continuum. 
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Fig. 3. — The spectral region around the [OIII] doublet for some of the objects included in 
the clean sample (see the main text). The spectra have been normalized with respect to the 
peak of the main [OIII] line. 
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Fig. 4. — Rest frame position of the two lines of the doublet of [OIII] in the clean sample 2 
(see main text). The vertical and horizontal lines indicate the theoretical local values. 
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Fig. 5. — Comparison between the estimations of Aa/a by each of the methods used in this 
paper (see main text for the details). 
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Fig. 6. — Aa/a vs. redshift. Each bin contains the contribution of ~ 175 spectra. 
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Fig. 7. — Histogram of the line ratio A of the [OIII] hnes. 
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Fig. 8. — Comparison between the estimations of Aa/a obtained from [OIII] and those 
derived from [SII] {left) and [Nelll] {right) respectively. 
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Fig. 9. — Rest frame position of the two lines of the doublet of [Nelll]. The vertical and 
horizontal lines indicate the theoretical local values. 
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Fig. 10. — Comparison between the values of Aa/a estimated by Bahcall et al. and in this 
study. Results from spectra included in the raw (clean) sample are plotted as empty (solid) 
points. 
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Table 1: Aa/a at different redshffts 



Redshift 


Aa/a (10-5) 


l-cr errors 


0.003 - 0.076 


+0.6 


5.2 


0.076 - 0.110 


-10.3 


5.7 


0.110 - 0.152 


+21.5 


8.0 


0.152 - 0.199 


+6.1 


7.4 


0.199 - 0.240 


+7.0 


6.5 


0.240 - 0.300 


-0.3 


8.6 


0.300 - 0.378 


-0.8 


8.4 


0.378 - 0.490 


+6.6 


7.5 


0.490 - 0.747 


+1.4 


8.2 



Table 2: Results for the different atomic species 

Doublet Spectra Aa/a (10~^) l-cr errors 

OIII 1568 +2A 2l 

NcIII 168 +358.5 11.2 

SII 481 -18.8 41.8 



